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Abstract: Consider a system of n polynomial equations and r polynomial inequations in 
n indeterminates of degree bounded by d with coefHcients in a polynomial ring of s parame- 
ters with rational coefHcients of bit-size at most a. From the real viewpoint, solving such a 
system often means describing some semi-algebraic sets in the parameter space over which 
the number of real solutions of the considered parametric system is constant. Following the 
works of Lazard and Rouillier, this can be done by the computation of a discriminant vari- 
ety. In this report we focus on the case where for a generic specialization of the parameters 
the system of equations generates a radical zero-dimensional ideal, which is usual in the 
applications. In this case, we provide a deterministic method computing the minimal dis- 
criminant variety reducing the problem to a problem of elimination. Moreover, we prove that 
the degree of the computed minimal discriminant variety is bounded by D :— {n + r)(i'^"+^' 
and that the complexity of our method is o-'^(i)d'^("+'*) bit-operations on a deterministic 
Turing machine. 
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Complexite de resolution d'un systeme parametre 
d'equations et d'inequations polynomiales 

Resume : On considere un systeme de n equations polynomiales et r inequations en n 

inconnues et .s parametres. Le degre des polynomes consideres est majore par d et leurs 
coefficients sont rationels de taille binaire au plus cr. D'un point de vue reel, resoudre un tel 
systeme revient souvent a decrire un semi-algebrique de I'espace des parametres au-dessus 
duquel Ic nombrc dc solutions reels du systeme parametre considere est constant. D'apres Ics 
travaux de Lazard et Rouillier, on pent obtenir ce semi-algebrique par le calcul d'une variete 
discriminante. Dans ce rapport, nous nous restreignons au cas ou le systemes d'equations 
donnes en entree est zero-dimensionel pour une specialisation generique des parametres, ce 
qui correspond a une situation courante dans les applications. Dans ce cas, nous proposons 
une methode deterministe pour calculer la variete discriminante minimale en reduisant le 
problcmc a un problcmc d'climination. Dc plus, nous prouvons que Ic degre de la variete 
discriminante minimale est majoree par D := {n + r)S"'^^^ et que la complexite de notre 
methode est de cr^^^^i)^^""'"*^ operations binaires sur une machine de Turing. 

Mots-cles : Systeme polynomial parametre, Variete discriminante. Elimination, Complexite 
deterministe 
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1 Introduction 

The parametric polynomial systems are used in many different fields such as robotics, op- 
timization, geometry problems, and so on. In (2^] the authors introduce the notion of 
discriminant variety which allows them to split the parameter space in open cells where the 
number of real solutions is constant . Even if it is efficient in a practical point of view, their 
algorithm is based Grobner bases computations, whose complexity is not yet well under- 
stood. Thus it does not allow us to give a better bound than the worst case's one, which is 
in exponential space (;19:). 

In this article we prove that, under some assumptions, the computation of the minimal 
discriminant variety of a parametric system is reducible to the FPSPACE problem of general 
elimination The proof of the reduction correctness presented here is non trivial. The 
reduction itself is simple and preserves the sparsity of the input system. 

Our input is a system of polynomial equations and inequations of degrees bounded d, 
which can be written as: 

fiit,x) = ( giit,x)^0 

: and i : {t, x) e C x C" 

fn{t,x)^0 [ gr{t,x)^0 

where x are the unknowns and t are the parameters. Moreover, for all specializations in an 
open ball of the parameters space, the system has a finite number of simple solutions in the 
unknowns. Such a system will be said generically simple (see Definitional- 
We prove that the degree of the minimal discriminant variety of a generically simple para- 
metric system is bounded by 

{n + r)d"+i 

Our algorithm for generically simple parametric systems runs in 

bit-operations on a deterministic Turing machine. 

When we aim to solve a parametric system, we face two kinds of issues: either we want to 
describe the solutions in terms of the parameters, or else we want to classify the parameters 
according to properties of the parametric system's solutions. Different methods have been 
developped to treat these two problems. 

Regarding the first one, many algorithms exist in the literature. Among them we may cite 
rational parametrizations ISD,, triangular sets decompositions (2^1, comprehensive Grobner 
bases |34l . We may also mention numerical algorithms such as the Newton- Raphson or 
the homotopy continuation method |32l I3H . which can be used after a specialization of the 
parameters. 
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Regarding the second problem on the parameters classification, few algorithm are avail- 
able, whereas many appHcations face it, such as parametric optimization (^3)' robot mod- 
eling (JU), geometry problems (^Q) or control theory ([2j) for example. The C.A.D. [Till 17] 
is the most widespread method. It computes an exhaustive classification, leading to a com- 
plexity doubly exponential in the number of unknowns. Some of the algorithms mentioned 
above ([33, 22j) may also return such kinds of classifications. Especially in j22] the authors 
compute a complete partition of the parameters space in constructible sets where the vector 
of multiplicities of the system's solutions is constant. The time complexity of their algorithm 
is d'-''-" However, they don't consider inequations and their algorithm is not meant to be 
implemented. The minimal discriminant variety is included in both of the precedent com- 
putations. It describes the maximal open subset of the parameters space where the system's 
solutions evolve regularly. The computation of this variety is indeed sufficient for a lot of 
appHcations. 

Our method is a reduction to the general elimination problem. The elimination problem 
has been widely analysed in the past decades, as it is a key step for quantifier elimination 
theory (in [23 EHl IH E] for example), computation of the dimension of an ideal (f6] among 
others) or impHcitization theory (see JL2j)- Different techniques and software have been de- 
veloped. We may mention sparse resultants (see ^31 and references therein), linear system 
reductions (in fSj for example) , linear systems parametrized with straight-line programs (see 
[28j 24J), parametric geometric resolution ( ^21. .30^ 1 or Grobner bases (see J2j and filil ilTi] 
for the last improvements). 

This article is divided in three parts. In the first one we reduce the problem of computing 
the minimal discriminant variety to the elimination problem. In the second part, we bound 
the degree of the minimal discriminant variety. And in the last part we give some examples. 

Definition and notation 

In the following, we assume that 



are some polynomials in degrees di — deg{fi) and d'j — deg{gj) for 1 < i < n and 1 < j < r. 
We denote by P„ the projective closure of C" and by tt : x C" ^ C (resp. 7f : C x P„ ^ 
C) the canonical projection onto the parameters space. The exponent h (resp. hi) of a 
polynomial or of an ideal denotes its homogenization by the variable Xq with respect to 
the variables Xi, - ■ ■ ,Xn (resp. its homogenization by the variable Xi with respect to the 
variables Xq, ■ ■ ■ , Xi, ■ ■ ■ , Xn) ■ The term parameters will refer to the variables Ti, ■ ■ ■ ,Ts, 
while the term unknowns will refer to the variables Xi, • • • , Xn- 

Finally we use the following notation for the specialization of some variable. For / C 
Q[yi, • • • , Yfc, Z] and a G Q, we denote: 



fir ■ ■ , fn, gi, 



,5, GQ[ri,---,r,][Xi,...,x„] 




—a 



{I+{Z~a))nQ[Yi,---,Yk] 
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In order to define the notion of discriminant variety according to our assumptions, we in- 
troduce the notion of geometric regularity. 

Definition 1 Let E be a subset of the parameters space. 

A parametric system S defining a constructible set C is said to be geometrically regular over 
E iff for all open setU C E, ir restricted to ir^^iJA) DC is an analytic covering. 

The minimal discriminant variety is now defined as follows. 

Definition 2 f2f)^ A discriminant variety of the parametric system S is a variety V in the 
parameters space such that S is geometrically regular over \ V . 

Among the discriminant varieties we define the minimal one: 

Definition 3 The minimal discriminant variety of S is the intersection of all the dis- 
criminant varieties of S. 

For the computation of the minimal discriminant variety, we will assume some properties 
on the input parametric systems we consider. 

Definition 4 Let S be the parametric system defined by: 

r fi{t,x)^0 ( gi{t,x)^0 

< ■■ and \ ■■ {t,x) e e X C" 

i fn{t,x)^0 [ 9r{t,x)^0 

Denoting 11^=1 93; assume that the ideal in the polynomial ring over the field of 

fractions of the parameters 

= (/i, •••,/«): 3^ C Q(Ti, ••• , T,)[Xi, • • • , Xn] 

is radical and 0-dimensional. 
Then S is said generically simple. 

Remark 1 Note that the ideal I generated by /i, • • • , /« C (Q[ri, • • • , T,, Xi, • • • , X„] needs 
neither to be radical nor equidimensional, although it is sufficient to satisfy the hypotheses. 

Moreover, given a parametric system S defined by /i = 0, • • • , /„ = 0, 31 7^ 0, • • • , .g^ 7^ 0, we 
introduce these two polynomials: 

— js is the determinant of the Jacobian matrix of /i, • • • , /« with respect to the un- 

knowns, of degree denoted by 6 

— gs is the product of the gi for 1 < i < r of degree denoted by S' 
Note that we have S < di — n and S' — J2'j=i ^'j- 
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Main results 

We can now state our main results. 

Theorem 1 Let S be a generically simple parametric system. 

Then the total degree of the minimal discriminant variety is bounded by 

di---dnil + S + S') 

Theorem 2 Let S be a parametric system generically simple defined by fi = 0, ■ ■ ■ , fn = 
0,gi 0, - ■ ■ ,gr ^ 0. Then the union of the varieties defined by the n + 2 following ideals: 
- R denotes the ring Q[Ti, ■ ■ ■ ,Ts] 

((A\ • • • , ZXog^s - 1, ^1 - 1) n R[Xo])^^^^^ (Xi) 
((/^ • • • , fH, ZXog'^s -l,Xn-l)n BIX,]) i^^^g (J„) 

((/li ■ ■ ■ T fn, 9s — Xn+1, ZXn+1 — 1) H ^[-'^n+l] ) |x„+i=0 i^n+l) 

{{fi,---Jn,js,Zgs~l))nR iIn+2) 
is the minimal discriminant variety of S . 

Corollary 1 A discriminant variety of a generically simple parametric system can be com- 
puted in: 

a'^«(di---d„(<5 + 5'))°("+^' 

steps on a classical Turing machine. The variable a denotes the maximal binary size of 
coefficients 0/ /i, • • • , /„ and 51, • • • , g^. 

Remark 2 If the system is not generically simple, then the the union of the varieties com- 
puted is the whole parameter space, which is thus an easy way to check if the initial conditions 
are verified. 

Remark 3 Any elimination algorithm may actually be used to compute a discriminant va- 
riety, which is welcomed when it comes to an effective computation. Among others, Grobner 
bases with a block ordering /J ,51 \16^ . sparse elimination ^ISJ or straight-line programs ]2f^ 
may lead to efficient computations. 

Remark 4 If we allow ourself to use the model of a probabilistic bounded Turing Ma- 
chine, then at the cost of the sparsity of the system, we may replace the computation of 
V(Xi), . . . , V(X„) by the computation of the variety of: 

{{ft • • • , - 1, 71^1 + ■ ■ • + inX,, - 1) n Q[ri, • • • , r,] [x,])\x„=o 

where (71, . . . , 7„) is chosen randomly in {0, 1}" and D Sdi • • • (i„.^ 

^The remark^ and the corollary Q are proved Section |3] 
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2 Log-space reduction 
2.1 Preliminaries 

The goal of this section is to show how to reduce the problem of computing the minimal 
discriminant variety (the discriminant problem) to the elimination problem. We know that 
the elimination problem is solvable in polynomial space (f^). Thus via the reduction 
we prove that the problem of computing the minimal discriminant variety is solvable in 
polynomial space. 



Discriminant Function: 
Input: fi, - ■ ■ , fn,gsjs G Q[Ti, - ■ ■ ,Ts,Xi, - ■ ■ ,Xn] 

Output. ^^^^ ^ ^ ■ ,qt,ut e Q[Ti, ■ ■■ ,Ts] such that U-^iV((gi4, • • • is the mini- 

mal discriminant variety. 



Elimination Function: 
Input: / e Q[Ti,---,T,][Xi,---,X„]; 

1 Ti,...,r, 

Output. . . . ^ ^ Q[Ti, • • • , T,] such that Y{{qi, ■■■ ,qt)) is the variety of the elim- 

ination ideal (pi, • ' ' ,Pm) nQ[Ti, ■ ■ ■ ,Ts]. 



To achieve the reduction, we will first describe more precisely how the minimal discrim- 
inant variety can be decomposed. In [221, the authors show that the minimal discriminant 
variety of a generically simple parametric system S is the union of 3 varieties, denoted re- 
spectively by Vinf , Vineq and Vcrit ■ Let us remind the definitions of these varieties under 
our assumptions. 

Definition 5 Let S be a generically simple parametric system defined by fi — 0, ■ ■ ■ , fn ~ 
and gi ^ 0, • • • , ^ 0. The varieties Vinf , Vineq and Vcrit of the parameters space are 
respectively defined as follow: 

Vnf =Tf{Cs_nnoo) 

where Cs is the projective closure of the constructible set defined by S, and Hoc = 
(C^ X P,i) \ (C^ X C") is the hypersurface at the infinity. 

Vneq = v{{is : 5^ + (gs>) n Q[ri, • • • , r,]) 
Vc„t = v{{is : 5^ + {is)) n Q[ri, • • • , t,]) 

Theorem 3 126] The minimal discriminant variety of a generically simple parametric sys- 
tem is the union of Vinf , Vineq and Vcrit ■ 
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Geometrically, this theorem characterizes the different varieties in the parameter space over 
which the generically simple parametric system is not geometrically regular. More precisely, 
the theorem means that over the minimal discriminant variety, three types of irregularity 
may appear. The first one is the intersections of the system of equations with the Jacobian. 
The second one is the intersection with the inequations. And the last one is the intersection 
in the projective space of the the hypersurface at the infinity with the projective closure of 
the parametric system's zeros. 

Vcrit is already directly the solution of an elimination problem. This is the component 
for which the generic radicality condition is needed. We will now focus on reducing the 
computation of each of the two varieties Vinf and Vmeq to the elimination problem. 

2.2 Reduction of Vinj and correctness 

Before going further, it should be clear that the computation of Vinf can not be handled 
by the standard projective ehmination methods if we want to certify a singly exponential 
complexity. All of these methods have no good complexity bounds essentially because of the 
intersection with the particular hypersurface at the infinity as we will see later. However 
this doesn't prevent us to use results of the projective elimination theory. 

Using the algebraic representation of the projection tF of [T^, with the notations of the 
definition we reformulate Vinf. 



where Js := {Is ■ Os')''- Note that Cs = V(Js). 

And using the reformulation of the ideal homogenization of [12], we obtain a formulation 
of Js which match expHcitly the input of the problem: 



This is however not yet satisfying since this formulation is not trivially reducible to a single 
elimination problem. The problem here does not come from the saturation by the variables 
Xi which can be simply handled with the Rabinowitsch trick (2^1 of adding the new variable 
Z and the new equation ZXi — 1 to the initial polynomials. Neither is the saturation by gs a 
problem since again we may add the equation Zgs— 1 = 0. The complications arise actually 
from the variable Xq. First we have to saturate by Xq and then we have to specialize Xq 
with to finally eliminate the variables Xi. And it is regrettable since this prevents us to 
use the usual trick to get rid of the saturation, as we saw in introduction. Moreover we 
don't want to apply successively two ELIMINATION Function since it could lead us to an 
exponential space algorithm. 

Fortunately we manage to sort out this problem by proving that for the variety we want 
to compute, we can commute the speciaHzation of Xq by and the elimination, which is 
remarkable since this operation will allow us to use the Rabinowitsch trick to localize by 




Js 







oo 
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Xq. Note that the commutation step does not alter the computation only because of the 
particular structure of Vinf ■ 

Proposition 1 Let S be a parametric system. Then the component Vmf of the minimal 
discriminant variety of S is the union of the varieties defined by the n following ideals for 
1 < i < n: 

{{ft ft ZX^g), - 1, X, - 1) n i?[Xo]) 

Remark 5 Note that the condition generically simple is not needed for the reduction of the 
computation of Vinf ■ Moreover the proposition remains true even if the number of equations 
differs from the number of unknowns. 

The proof of this proposition is based on the three following lemmas. The first one gives 
some basic useful equalities, where hi denotes the homogenization by the variable Xi with 
respect to the variables Xq, • • • , X^, • • • , X„. 

Lemma 1 Let J C Q[ri, • • • , Ts] [Xq, • • • , X„] be an ideal homogeneous in X^, • • • , Xn 
and p be a polynomial of Q[Ti, ■ ■ ■ , Ts] [Xq, • • • , Xn] also homogeneous in Xq, • • • , Xn- Then 
for all < i < n we have: 

(Jix.=i)'- ^j--xr 

{J : pn\x^=i ^ J\x^=i ■■ pTx.=i 
J : n Q[ri, • • • , Ts\ = J|x,=i n Q[Ti, • • • , T,] 

and for all 1 < i < n: 

n Q[Ti, • • • , Ts][Xo] = (J n Q[Ti, • • • , Ts][Xo,X,])\x,=, 
Proof: These are classical results that can be recovered from 12j. □ 

Now comes the first lemma toward the reduction, which proves essentially that the union of 
the varieties defined by the ehmination ideals of the proposition contains Vinf- 

Lemma 2 Let J be an ideal o/Q[Ti, • • • , Tg] [Xq, • • • , Xn] homogeneous in Xq, • • • , Xn- Then, 
for all 1 < i < n we have: 

(jnQ[ri,...,T,][Xo,x,])|Xo=o,x,=i 

n 

(J|Xo=o :^r)nQ[Ti,...,T,] 

Proof: Let p e (>/ H Q[ri, • • • , Ts] [Xq, Xi])|Xo=o,Xi=i- The polynomial p is homogeneous in 
Xq, ■ ■ ■ , Xn since it depends only on the variables Ti, . . . ,Ts. Thus with the notations of the 
lemmaQl we have p € {{J\Xo=o) 1X1=1)'''' ■ And J|Xo=o being homogeneous in Xq, ■ ■ ■ ,Xn, 
one can apply the first equality of Lemma to deduce p G J\Xo=o ■ which proves the 
desired result. □ 

And finally comes the keystone lemma related to the proposition, proving the reciprocal 
inclusion. 
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Lemma 3 Let J be an ideal o/Q[Ti, • • • , Ts][Xa, • • • , Xn] homogeneous in Xq, • • • , X„. Then, 
for all 1 < i < n, we have: 

^(jnQ[ri,...,T,j[Xo,x,j)|x„=o,x,=i 
u 

n;=i(^|Xo=o:^r)nQ[Ti,...,T,] 

Proof: Let p G r\]=iiJ\Xo=o ■ X°°) n Q[Ti, • • • , T^]. By definition there exist 91, • • ' , Qn G 
Q[Ti, ■■■,Ts] [Xo, • • • , X„] and fci, • • • , A:„ G N such that: 

' pi := pX^^ + Xoqi 

: G J 

^ Pn pX^" + Xog„ 

Since the part of Pi of degree h in Xo,---,X„ belongs also to J, we can assume that 
pi, - ■ ■ ,Pn are homogeneous in Xq, • • • , X„. Thus, we have in particular: 

degjfi,...^x„fe) < kj 

Now we fix a total degree term order <x on the variables Xi, - ■ ■ ,X„. Let K denote the 
field Q{Ti, ■ ■ ■ ,Ts, Xq) and consider Pi, - ■ ■ ,Pn as polynomials of K[Xi, ■ ■ ■ , Xn]. Denoting 
by J the ideal they generate, it follows immediately that 

G := {pi, ■ ■ ■ ,Pn} 

form a Grobner basis of with respect to <x since the pi have disjoint head terms. Let 
i be an integer between 1 and n. We first show how to prove the lemma when we have a 
polynomial of J such that: 

- it is univariate in (1) 

- it has all its coefficients in Q[ri, • • • , T^, Xq] (2) 

- its head coefficient is a power of p (3) 
Assume for a while that is such a polynomial, dx^ being its degree in Xi. It follows indeed 
that ri G J''^ the contraction ideal of J^. And since p = lcm{HC{g)\g G G} we have 0: 

J'' - (G> : 

meaning that for some k £ N, p'^Vi G (G) C J. Finally J is homogeneous so that f^, the 
part of degree dx^ of p^Vi, belongs to J n Q[Ti, • • • , T^] [ATo, Xi] and can be written as: 

n = p'xf^ + x^q 

with I G N and q G Q[Ti, • • • , Ts][Aro, Xi], which is an equivalent way of writing 

pe^iJn Q[Ti, • • • , T,][Xo, X,])\x,=o,x,=i 
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It remains us to show the existence of a polynomial satisfying (1),(2) and (3). To carry out 
this problem, we first notice that J is zero-dimensional in K[Xi, ■ ■ ■ ,X„] since the set of 
the head terms of its Grobner basis contains a pure power of each variable Xi. So, we may 
consider the finite dimensional /iT— space vector A = K[Xi, ■ ■ ■ ,Xn]/J along with e the 
monomial basis of A induced by G. More precisely, denoting by x the class of x in A, we 
define see e as the set of ey for 1 < j < I? := d\m{A) such that Cj is a term of K[Xi, • • • , X„] 
not multiple of any head term of G. Finally we denote by S the multiplicatively closed set 
{p'^,k G N}. We will follow a classical method to exhibit a monic univariate polynomial 
from a zero-dimensional ideal, with coefficients in K. And with results of we ensure that 
its coefficients are not only in K but rather in the ring S~^Q[Ti, • • • , Ts,Xo] C K. 
Let us introduce the classical linear application of multiplication by Xi : 

<Pi: A ^ _A_ 

q ^ XiQ 



Then we note M, the matrix of in the base e: 



M, 




we notice that the coefficients of Mi come from the reduction of the monomials XiCi for 
1 < / < Z) by the Grobner basis G. And as we can see in j^], this kind of reduction only 
involves division by the head coefficients of G, such that: 



Xiei = ci^ieij H h CD,ieD,i 

with ci^i, • • • , cd,i not only in K but more precisely in the ring 5^^Q[ri, • • • , Ts,Xo] C K 
where S = {p^,k e N}. As a straightforward consequence, if we denote by Vi the monic 
characteristic polynomial of in the new variable ?7, we have Vi € S'~^Q[ri, • • • , T^, Xo][C/]. 
Besides by the Cayley-Hamilton's theorem, Vi applied to the variable Xi is the null element 
of A, meaning that Vi{Xi) belongs to J and may be written as: 

V^{X,) = XP + Cd-iX^-' + --- + Co 
with Cfc e S'^^QlTi, • • • , T,, Xo] for 1 < fc < £> - 1. Finally, for some fc' e N we have 

r, ■.^p'''v^iX^) ejnQ[Ti,---,T,,Xo][X,] 
which satisfies all the conditions we wanted to achieve the demonstration. □ 
Finally, a proper combination of the lemmas proves the proposition 
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2.3 Reduction of Vineq and correctness 

As to bound the computation of the variety induced by the inequations 

Vine, = V ((7s : gf + (gs)) n Q[Ti, • • • , T,]) 

the direct approach consists first in performing a saturation and then in using the output 
along with gs as the input of an elimination algorithm. However this method may not have 
a single exponential bound on the time complexity in the worst case. Hence both of these 
algorithms may use a polynomial space in the size of the input, which could finally cost an 
exponential space if no more care is taken. 

In this section we show how to bypass the problem, notably by relaxing the condition on 
the output and allowing some components of Vinf to mix in. 

Proposition 2 Let S be a parametric system. If we denote by Wineq C the variety 
defined by the following ideal: 

((/l: ■ ■ ■ J fn, gs — ^n+l, ZXn+i — 1) 

nQ[Ti,...,T«][X„+i])|^„^^=o 
then the following inclusions chain holds: 

Vineq C W^ineq C Vineq ^ Vinf 

The first step to prove this proposition is to delay the saturation. 

Lemma 4 Let pi , • • • , p„i, q,r G Q[Yi , • • • , 11] . Let us fix < a term order and assume that 
the head monomial of q shares no variable in common with the monomials of pi, ■ ■ ■ jPrntf. 
Then we have the following equality: 

{Pi, ■ ■■,Prn) ■■ r°° + (q) = {pi, ■ ■ ■ ,Pm,q) ■■ r°° 

Proof: The inclusion from left to right is trivial. For the other inclusion, let p £ (pi, • " " iPmi q) '■ 
r°°. Denoting by M the head monomial of q with respect to <, we obtain by division: 

p = p' + qt p',t€q[Yi,---,Yk] (1) 

such that no monomial of p' is multiple of M. It remains to show that p' belongs to 
(pi, • • • ,Pm) '■ and the proof will be complete. By hypothesis, we know that there exists 
I > and Ci, • • • , Cm, c G Q[Yi,- ■■ ,Yk\ such that: 

r'p' = cipi H h CmPm + cq 

We divide each of the Cj by q as in (1) and denote by the remainder of the division. We 
thus obtain: 

rV - c'lPi c'^pm = c'q 
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with c' £ Q[Yi, • • • , Yfe]. We remark that the polynomial on the left part of the equality 
has no monomial which is multiple of M, while the head monomial of the right part of the 
equality is M times the head monomial of c', which means c' = and this achieves the 
proof. □ 



Corollary 2 Let fi, ■ ■ ■ , fn, 9 be some polynomials of Q[Ti, ■ ■ ■ , Tg] [Xi, • • • , X„] . Then: 

ifl, • • • , /n) : 5°° + (5 - ^n+l) - (/l, • • • , /n, 5 - ^n+l) : ^^+1 

Thanks to this result, we can now reformulate Vineq as being the variety of: 
((/i, ■■■Jn,gs- ^„+i) : ^r+i) |x„^,=o n 'Qt^i' • • • ' 

The reduction is not yet complete and we encounter here the same problem we had for the 
computation of Vinf, that is the saturation by Xn+i before the specialization of Xn+i by 
0. This is just fine since the lemmas fll2l and l3l provide us tools to handle it, even if they do 
not completely solve the problem yet. 
For the first inclusion, we note: 

■— ifl, ■ ■ ■ , fn, gs — Xn+i) : 
it follows that the varieties of the proposition [2 rewrite as: 

w,„e, - v((/'5nQ[ri,...,r,][x„+i])|x„^^=o) 

and we show easily: 

(/^ n Q[ri, • • • , r,][x„+i])|x„^^^o c /|i„^,^o n Q[Ti, • • • , t,] 

which, in term of varieties, proves the first inclusion of the proposition [2 

For the other inclusion we will mainly use the lemma El For this, we introduce the 
homogenization variable Xq, and denote with the exponent h the homogenization by Xq 
with respect to Xq, • • • , Xn+i- We need also the following classical lemma, which dissociates 
the affine part from the component at the infinity of a homogeneous ideal. 

Lemma 5 Let J G Q[Ti, - ■ ■ , Ts] [Xq, ■ ■ ■ , Xn+i] an ideal homogeneous in Xq, ■ ■ ■ , Xn+i- 
Then the following equality holds: 

Vj^VJ+ i^o) n y/J : x^ 
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In term of varieties, the equality follows from the observation that V( J) is the union of 

V( J) n Hoo and V(J)\Hoo. 

We now homogenize by Xq and we get: 

Using lemma we get directly: 

Then we show that (/'^'')[x„+i=o contains an ideal which begins to look like what we want: 

I'' - :gr + (.9s-^„+i))" 

D Js + (Xo.9§) 
Then, the lemma fallows us to split the ideal Js + (^o.9s) i^- 



voc 
^0 



h h 
such that we now have the following inclusion: 

Wi„eg C V I fl (/i : Xf) n [h : ) n Q[Ti, • • • , T,] 

From there, denoting again (Q)[Ti, . . . , T^] by R, we remark for < J < "-: 

h■.X°°f^R D (Js)|jf„=o ni? 

And: 

/2 : n i? D (J5 + (.9s)) : : Xo°° n i? 

3 ((J5)|Xo=i + (.9s)) iXfni? 

D l'^ ■.X'^C]R 

Which allows us to conclude with: 

w„,e, c V (/^^ n n;=o (^s|xo=o) : xf) n • • • , r,] 

This proves the theorem [2 
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3 Degree issues 

The study of the degree of the minimal discriminant variety reUes strongly on the Bezout- 
Inequality jlHlEl- What we call degree of an ideal / (resp. a variety V) and denote deg(J) 
(resp. deg{V)) is the sum of the degrees of the prime ideals associated to -s/T (resp. the sum 
of the degrees of the irreducible components of V). With this definition, from |2I^I 118) we 
have for /, J C Q[Ti, . . . ,T„Xo, . . . , X„] and / G Q[Ti, . . . ,T„Xo, . . . ,X„]: 

deg(/ + J) < deg(/)deg(J) 
deg(/:/°°) < deg(/) 

deg(/nQ[ri,...,r,]) < deg(/) 

deg(/) = deg(V(/)) 

Degree of Vinj 

Here we use the prime decomposition of \/Js to bound the degree of Vinf ■ This decompo- 
sition will also allow us to prove Remark 31 
First we remind that from proposition 




fl Js :xr nQ[ri,---,r„Xo] 



where Js-(/^•••,/^) : 

Continuing with the properties of the degree we have: 



''OO 



deg(J5)<deg((/^,...,/^)) <di...d„ 
Let ^i, . . . be the prime ideals associated to \/Js. Then we have: 

<Pin •••n'Pfc = yjs 

deg(^Pi) + --- + deg(<Pfc)<di---d„ 

Now let denote by Ai, . . . , Aj the indices of the prime ideal which do not contain any power 
of Xi for some 1 < i < n. It follows that: 

n 
i=l 

such that 

deg{V,nf) = deg (niLi ^/J~rxf^Q[Ti, • • • , T., Mxo=o) 

<di---dn 
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We use the decomposition of yOs to prove the remark 31 
Proof: (of Remark^l 

We extend Lemma where we replace Xi by a homogeneous linear form in Xq, ■ ■ ■ 
which leads to the following property. If J is an ideal of Q[Ti, • • • , Ts][Xo, • • • , X„] homoge- 
neous in Xq, • • • , X„, and L € Q[^oj ' ' ■ > ^n] is a homogeneous linear form in Xq, ■ ■ ■ , X^, 
then: 

J:L°°nQ[Ti,...,T„Xo] = [J + {L - 1)) n Q[T,, ■ ■ ■ ,T,, Xq] 

From there, we know that the prime ideals which contain a power of Xi for all 1 < i < n 
contain in fact all the homogeneous linear forms of Q[A'o, • • • ,Xn]. Let denote by E the 
Q-spacevector of homogeneous linear forms of Q[^o, • • • , Xn]- Thus we have for all L e E: 

j 

i=l 

Let B denote the bounded lattice {0, . . . ,D — 1}" of E, where D — Sdi ■■■<!„. And A be 
defined by: 

A:=|J(<P;,,ni?) 

i=l 

Such that for L e B \ A, we have: 

And since each n i? is a strict linear subspace of E, it follows that A is the union of 
j — n"=i '^i = T strict linear subspaces of E . Each ^Pa; H E intersects the lattice B in at 
most Z?""^ points. Thus the probability of choosing L in B (1 A is < |. And for all 

L ^ B \ A we have: 



v^nf = V f|*p,;:L°°nQ[ri,-..,r„Xo]|x„=o 



' / k 



V 'P, + (L - 1) j Q[ri, • • • , T„Xo]|Xo=o 

V(((/^,...,/;^,ZXog§-l,L-l) 

nQ[ri,-.-,T„Xo])|Xo=o) 



□ 
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Thus with the properties of the degree, we have respectively: 



deg{Vineq) < dl ■ ■ ■ dn6' 

deg{Vcrif) < di---dn5 



Hence we proved the theorem 

Degree of representation of the elimination 

To compute the ELIMINATION Function in a deterministic way, we follow the ideas of 
which uses the afEne effective Nullstellensatz to reduce the problem to a linear algebra 
system of non homogeneous linear form. One could use the ideas of [28, 20, 21J to perform 
this elimination, whose complexity bounds rely on the projective effective Nullstellensatz of 
|25] . However these bounds only hold in a bounded probabiHstic Turing machine. 

Here we will use the Brownawell's prime power version of Nullstellensatz (see 0), which 
is a variant of the afHne effective Nullstellensatz: 

Theorem 4 |^ Let J C k[xQ, • • • , Xn] be an ideal generated by m homogeneous polynomial 
of respective degrees d2 > ■ ■ ■ > dm > di and 971 — {xq, • • • , a;„). Then there are prime ideal 
?Pi, • • • , *Pr containing J and positive integers eo, • • • , Cr such that: 



where /i = min(TO, n) 

Using the proposition 3 of j22], we know that if is a prime ideal, then there is n + 1 
polynomials /i, ■ • • , fn+i such that: 



r 



eo + Y.e,degi%) < (3/2)^di • • • 



i=l 



V(/i, 



with deg(/i) < deg(q3) for all 1 < i < n 
Thus we deduce the following: 



+ 



1 
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Proposition 3 Let I C Q[Ti, ■ ■ ■ ,Ts] [Xi, • • • , Xn] generated by fi, ■ ■ ■ , fm indexed such that 
their degrees satisfy d2 > ■ ■ ■ > dm > c?i- Then, with /i = min(m, n) we introduce: 

( m g,eQ[Ti,---,T,][Xi,---,Xn 
F := < ^gjtl and 

deg{g^f^)<{3/2)^'dl■■■df, 

Then we have: 

v(/ n Q[Ti, • • • , r,]) = v(^^ n q[Ti, • • • , t,]) 

Proof: We homogenize the polynomials /i, • • • , /m by with respect to Ti, • • • , T^, Xi, • • • , Xn, 
and denote by J the ideal they generate. Then with • • • , being prime ideals containing 
J and verifying the theorem of Brownawell, it follows that the result holds when intersect- 
ing J and • • • , *Pr by Q[Ti, • • • , Tg, H]. Finally we use the Heintz's proposition reminded 
above on each *Pi and specialize by 1 to conclude. □ 

Now consider the coefficients of the polynomials gi,..,gm,g as unknowns. Assume further- 
more that gi, - ■ ■ ,gm contains all the monomials in Ti, ■ ■ ■ ,Ts, Xi, ■ ■ ■ , X„ of degree less or 
equal to (3/2)^di • • • d^, and that g contains the monomials in Ti, • • • , only. Thus, finding 
the coefficients satisfying the formula: 

m 
i=l 

reduces to the problem of finding null space generators of a matrix of size lower or equal to 

(m + l)((3/2)^di • • • X {{3/2rdi ■ ■ ■ 

Hence the complexity of the corollary follows. 



4 Example 

We show here an example of minimal discriminant varieties application in our framework. 
It will allow us to prove that the real parametrization of the Enneper surface matches its 
real implicit form. In the author solves this problem with a combination REDLOG,QEPCAD 
and QERRC. Through the process, he has to simplify formulas whose textual representation 
contains approximatively 500 000 characters. We will see that our framework allows us to 
use minimal discriminant varieties to solve this problem. Notably, this allows us to keep 
formulas small. The following computations are done with the Maple package DV, which 
uses FGb to carry out the elimination function. We also use the factorization functions of 
Maple to take the square-free part of the polynomials given in the input, and to simplify the 
output. Finally RS and the Maple package RAG allows us to treat the discriminant varieties 
we compute. All these software are available in the Salsa Software Suite p. 
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When E and F are two lists of polynomials, T a list of parameters and X a list of 
unknowns, we denote by 

T)Y{E,F,T,X) 

the discriminant variety of the parametric system S : (p = 0)p^E A (q 7^ 0)q^F- 



4.1 Definition of the Enneper surface 

The real Enneper surface £ C M'^ has a parametric definition: 

£ ^ {{x{u,v),y{u,v),z{u,v)) \ (u,w)eR^} 

x{u,v) = 3u + 3uv^ — u'^ 
y{u, v) = 3v + 3u^v — 
z{u, v) = — 

We will also consider the graph of the Enneper surface £g C defined as follows: 

£g ^ {{x{u,v),y{u,v),z{u,v),u,v) \ (m, u)eK^} 
Beside, a Grobner basis computation returns easily its implicit Zarisky closure £ |12lll4j : 

£ = {{x,y,z) eR^ I p{x,y,z) = 0} 

p{x, y, z) = -19683x6 + 59049x^^2 _ logasx^z^ - 118098x'*z2 + 59049a;'*z - 59049x^2/4 
-56862x^^2^3 „ ii8098x2y2^ - 1296x2^6 - 34992x22^ - 174960x22-^ 
+314928x^23 + 19683?/S - imiby^z^ + 118098?/'^z2 + 59049?/''z + l2%Qy'^z^ 
-34992^2^5 + 1749602/2^4 + 314928^^^3 _^ q^^^'s _ i0368z^ + 4199042^ 



4.2 Discriminant varieties 

The main idea to compare and £ is in a first step to compute the union of their discrim- 
inant varieties, V . In a second step we compare £g and £ on a finite number of well chosen 
test points outside of V . Finally, the properties of the discriminant variety ensure us that 
the result of our comparison on these test points holds for every points outside of V . 

More precisely, £g and £ are both algebraic varieties of dimension 2. Thus we choose a 
common subset of 2 variables, x and y for example, which will be the parameters for the 
two discriminant varieties: 

V^"^ := DV( [x - x(u, v),y- y{u, v), z - z{u, «)],[], [x, y] , [z, u, v] ) 



V^y-.^J^Yi [p(x,2/,z)] A]A^.V] 



\z\ 



The number of equations equals the number of unknowns in both case and our algorithm 
returns a non trivial variety for both systems. This ensures us that the two systems are 
generically simple. Here are the results of the computations, which lasted less than 1 second 
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on a 2.8 GHz Intel Pentium cpu: 

V^'^ = V(?/6 + QOy^ + 768y2 _ 4096 + 3x^y^ - 312x^y^ + 768x^ + 3x^y^ + 60x^ + x^) 
UV(a;6 + A8x^ + 3x^y^ - 336x^y^ + ix^y^ + 768a;2 + 4096 + 768y2 + 482/"* + y^) 

V"^ = V(?/6 + my^ + 768y2 - 4096 + 3x^y^ - il2x^y^ + 768x2 _^ ^^iy2 _,_ qq^a ^6) 
UV(x6 + 48a;'* + Scc^'y^ _ 336x^^2 + Zx^y^ + 768^2 + 4096 + 768j/2 + 48^-* + y^) 
UV(x - y) U V(y) U V(x + y) U V(x) 

We denote by -Kxy : M? ^ M? the canonical projection. Then the properties of the 
discriminant variety ensure us that for each connected component C of \ {y^v y V2^), 
{TT~y(C) n £, TTxy) and {'n~y(C) n 5, iTxy) are both analytic covering. Moreover, E <Z E. Thus 
if C is a connected component of M? \ [V^^ U ¥2^), we get the following property: 

3p e C, n £ = nS^ypeC, n^^ip) n f = ^-^^(p) n £ 

This allows us to prove that £ and £ are equal above ^ (y;^v y y^^^); we take one 
point p in each connected component of \^ (y^^ y Vj^*'), and check that the number of 
real solutions of T^^y{p) H £ and of n^yiji) n £ is the same. We use the RAG package to get 
one point in each connected component and RS to solve the corresponding zero dimensional 
real systems. This allows us to prove that 

n-y\R' \ U V2^y)) n£^ 7r-y\R^ \ (Vf ^ U y^"')) n £ 

In order to get more information, we repeat this process using respectively the dis- 
criminant varieties on the parameter set {x,z} and {y,z}. This leads to the following 
computations: 



and 



yxz 


:= DV( [x - x{u 


v),y-y{u,v),z~ z{u,v)] , [ 


] , [x, z] 


[y,u,v] 


) 




:== DV( 


[p{x,y,z)] , [ 


] , [x,z] 


[y] 


) 


vr 


:= DV( [x - x{u, 


v),y-y{u,v),z- z{u,v)] , [ 


] : [y, A 


[x,u,v] 


) 


yyz 


:= DV( 


[p{x,y,z)] , [ 


] , [y,z] 


[x] 


) 



The result is shown on Figure Q 

Then we compute as above one point in each connected component of the complementary, 
and this allows us to prove that: 

7r,'-,i(R2 \ (vr u v^n) n £ = n-H^' \ [v^ u v^D) nf 

and 

71^,1 (r2 \ (vr uvy'))n£ = Ti-} (r2 \ {vr u vD) n £ 

Using the following notations: 
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it remains us to check what happens above each component of 

yxy p yxz p yyz 

An idea is to set apart the linear components from the others. We introduce 
Vl ■■= V{x + y)U V(x -y)U V{x) U V{y) U V(z) 

and denote respectively V^^ \ Vl,V'-''' \ Vl and \ by V^, and V^. Using the 
RAGlib, we verify that 

V'-^y n v^^ n vy^ 

has actually no real points. 

It remains us to check what happens on each of the 5 linear components of Vl- The 
intersection of £g or £ with a Hnear component P may be seen as a Hnear substitution of a 
variable. This operation produces 5 pairs of varieties of dimension 2 (TableOJ. To check their 
equality, we use the same strategy as above and compute the 5 discriminant varieties with 
1 parameter, 3 unknowns of Ki^ . . . , if 5, respectively Vki , ■ ■ ■ , Vk^, , and the 5 discriminant 
varieties with 1 parameter,! unknown of Li, . . . , L5, respectively Vl^ , Vl^ . We check 
that Ki = Li for each point by connected component of the complementary of Vk^ U Vl^ , 
in less than 1 second. And at last we intersect again the varieties with their discriminant 
varieties, which reduces the problem to compare 5 pairs of zero dimensional systems. Thus 
we check that the equality holds for the finitely many points considered. Finally this allows 
us to conclude that £ = £. 

5 Conclusion 

We provided a deterministic single exponential bit-complexity bound for the computation 
of the minimal discriminant variety of a generically simple parametric system. Note that 
the complexity of our algorithm relies on the elimination problem's complexity. Thus in a 
probabilistic bounded Turing machine, the work of for example leads to a polynomial 
complexity bound in the size of the output. Or if we are only interested in the real so- 
lutions, then the use of the single block elimination routine of [HE] improves directly the 
deterministic complexity bound of our method. 

The reduction presented in this article is easy to implement in conjunction with a software 
performing elimination, as those used in ^^E], or j2l] for example. 

It would be worth studying the complexity of the computation of the minimal discrimi- 
nant variety when we have more equations than unknowns. 
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^9 


w 


Y{x) 


V(t/) 


V{z) 


V{y + x) 


V{y - x) 




Ki 


K2 




Ki 


K5 




f — v) 


( X — x{u, v) 


( X — x{u, v) 


( X — x{u, v) 


( X — x{u, v) 


System 


< y~y{u,v) 


- y{u,v) 


I y-y{u,v) 


l-x- y{u,v) 


\ X y(u, v) 




[ z — z{u, v) 


[ z — z{u, v) 


[ — z{u, v) 


(z — z{u, v) 


(z — z{u, v) 


Parameter 


z 


z 


X 


x 


X 


Unknowns 


y,u,v 


X, U, V 


y,u,v 


z, u, V 


Z, U, V 


Minimal 


Vk, = 


Vk, = 


Vks = 


Vk, - Vk, = 


Discriminant 
Variety 


V{z)UV{z-3) 
uy{z - 9) 


V{z)UV{z + 3) 
uy{z + 9) 


V{x) U V(a;2 + 2) 


V{x + 4) U V(x - 4) U V(a;2 - 8) 
UV(a;2 + 2) 





£ 


w 


Vix) 


V(y) 


V(z) 


V(y + x) 


V(y - x) 


saw 


Li 


L2 




L4 




System 












{sq.fr = 


sqfr{p{0,y, z)) 


sqfr{p{x,0,z)) 


sqfr{p{x, y, 0)) 


sqfr{p{x,-x,z)) 


sqfr{p{x,x,z)) 


squarefree) 












Parameter 


z 


z 


X 


X 


X 


U nknown 


y 


X 


y 


z 


z 


Minimal 


Vl, = y{z + 9) 


Vl, = Y{z - 9) 


Vl3 = 


Vl, = Vl, = 


Discriminant 


UV(z) U V{z - 3) 


UV(z) U Y{z - 3) 


V(x) 


V{x + 4) U V(a; - 


- 4) U V(a;2 - 8) 


Variety 


U\{z - 9) 


UV(^; + 9) 




UV(a;) 



Table 1: Discriminant varieties of the sub varieties 
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